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Abstract 



in 
o 
o 

The decays v e —>■ e~xf, i>* e — * e + X\ , when kinematically allowed, constitute 
a source of 100% polarised charginos. We study the process e + e~ — ► v* e v e — > 
<^ \ e + Xi£~xfi with subsequent chargino decays Xi ~~ > ^/xA^Xij Xi QQ'Xi or 

CN \ their charge conjugate X\ ~^ QQ'Xii Xi ~^ The kinematics of this 

process allows the reconstruction of the sneutrino and chargino rest frames, and 
thus a clean analysis of the angular distributions of the visible final state products 
in these reference systems. Furthermore, a triple product CP asymmetry can be 
, built in X\ hadronic decays, involving the chargino spins and the momenta of 

the quark and antiquark. With a good c tagging efficiency, this CP asymmetry is 
, quite sensitive to the phase of the bino mass term M\ and can test CP violation 

o 

a: 

Qh' 

& 1 Introduction 



in the neutralino sector. 



An international linear e + e collider (ILC) offers an ideal environment for the study 



of supersymmetry (SUSY) [1,2], if this theory is realised in nature. With adequate 
choices for the centre of mass (CM) energy and beam polarisations, the various pro- 
duction processes simultaneously present in e + e~ annihilation can be disentangled. It 
is then possible to measure sparticle masses, couplings and mixings [3], allowing for the 
determination of the Lagrangian parameters. The quantum numbers of the new parti- 
cles must be investigated as well, in order to confirm that they are the superpartners 
of the Standard Model (SM) fields. In particular, it is crucial to determine the spins of 
the SUSY particles. Spin-dependent effects serve not only to deduce the particle spins 
but may also be used to verify the predictions of the theory, using the input from other 
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experiments. It is also conceivable that spin distributions and asymmetries, if precisely 
measured, can be used to improve the measurement of some of the mixing parameters 
of the Lagrangian. 

The existence of purely supersymmetric CP violation sources is other of the aspects 
that require a thorough investigation. At present, all CP violation observed in K and 
B oscillations and decay (with the possible exception of the time-dependent CP asym- 
metry in — > 4>Ks, whose experimental situation is still unclear) can be explained by 
a single CP-violating phase in the Cabibbo-Kobayashi-Maskawa mixing matrix. How- 
ever, the CP breaking induced by this phase cannot account for the observed baryon 
asymmetry in the universe. The most general Lagrangian of the minimal supersymmet- 
ric Standard Model (MSSM) contains a large number of CP violation sources, which 
can potentially yield observable effects at low and high energies. Within the neutralino 
sector, CP-violating phases can appear in the bino and Higgsino mass terms, M 1 and 
/i, respectively. At low energies the phases 0i, M of these two parameters generally 
lead to unacceptably large SUSY contributions to electric dipole moments (EDMs), if 
they are different from 0,7r. Present constraints from EDMs place strong restrictions 
on the values of and but do not necessarily require that Mi and \i are real, and 
cancellations among the different SUSY contributions can occur [4-6]. This possibility, 
and the need for further CP violation sources beyond the SM, makes the investigation 
of CP breaking in the neutralino sector a compelling task to be carried out at a linear 
collider. 

In this paper we focus on the determination of the spin and spin-related properties of 
sneutrinos and charginos, including CP-violating spin asymmetries in chargino decays. 
We study sneutrino pair production in e + e~ annihilation at a CM energy of 800 GeV, 
as proposed for an ILC upgrade. We concentrate in the channels 

e + e~ -> u*u e -> e + Xi e'xt ~> e+ ^~X°i e~qq'x°i , 

e + e~ -> u*u e -> e + Xi e~xt -»• e + qq'x°i e~v^ + x\ , (1) 

with q — u,c, q' — d, s. This process has three advantages for our purposes: (i) the 
charginos produced are 100% polarised, having positive helicity in the decay z>* — > e + Xi 
and negative helicity in u e — > e~xt', (^) the kinematics allows the determination of the 
sneutrino and chargino rest frames, and then the study of angular distributions in these 
reference systems; (Hi) its cross section is large, and backgrounds with 5 energetic final 
state particles plus large missing energy and momentum are small. Muon sneutrinos are 
also a source of polarised charginos, but the cross section for z>*z> M production is much 
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smaller. Beam polarisations P e - = —0.8, P e + = 0.6 can be used to increase the signal 
cross sections but they do not affect the angular distributions and CP asymmetries 
studied. 

For definiteness, we consider a SUSY scenario like SPSla in Ref. [7] but with a 
heavier slepton spectrum (and nonzero phases 0i, so that sneutrino decays to 
charginos are kinematically allowed. The analysis of e~ , e + angular distributions in 
the z> e , v* rest frames provides a strong indication that sneutrinos are scalar particles 
and charginos have spin 1/2. The angular distributions of the Xi (xt) decay products 
i <?7 <?' Q.i <f) with respect to the x± (xf) s P m can be precisely measured as 
well. Furthermore, summing Xi an d Xi hadronic decays, a CP asymmetry based on 
the triple product s± ■ (p^ x p q2 ) can be built, where s± is the spin direction of the Xi 
or xt, an d p q2 are the three-momenta of the antiquark and quark resulting from its 
decay, respectively, in the xt res t frame. A good c tagging efficiency, as expected for a 
future linear collider, is essential for the determination of quark angular distributions 
and the CP asymmetry. For the latter it is necessary to distinguish between the quark 
and antiquark produced in the decay. This can be done combining c tagging and the 
knowledge of the final state muon charge. 

We note that in chargino pair production e + e~ — > XiXi the charginos are po- 
larised as well [8], allowing for the measurement of decay angular distributions and 
CP-violating asymmetries. However, in this process the momenta of the decaying 
charginos cannot be determined with kinematical constraints, and thus the analysis 
of angular distributions is less clean. Moreover, the backgrounds for XiXi m the 
semileptonic decay channel are huge, being e + e~ — > W + W~ — > d^vjj, I = e, fi the 
most troublesome, with a cross section (including ISR and beamstrahlung corrections) 
of approximately 3.5 pb at CM energies of 500 and 800 GeV, for P e + = 0.6, P e - = —0.8. 
At any rate, in SUSY scenarios where sneutrino decays to charginos are not kinemat- 
ically allowed, XiXi production seems the best place for the study of chargino spin 
properties and CP asymmetries in chargino decays. 

This paper is organised as follows. In section 01 we discuss sneutrino production 
and decay, focusing on the features most relevant for our work, and fix the SUSY 
scenario used. The procedure for the Monte Carlo calculation of the processes and 
reconstruction of the signals is outlined in section El In section 0] we introduce the 
different distributions analysed and present our numerical results. In section El we 
compare with results obtained from other processes and draw our conclusions. 
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2 Production and decay of sneutrino pairs 



Sneutrino pairs are produced in e + e~ collisions through the Feynman diagrams depicted 
in Fig. [TJ Their two-body decays v e — > e~xt, K ~~ * e+ Xi are mediated by the v e e Xi 
vertex which, neglecting the electron Yukawa coupling, is given by 

£» e ext = ~9Vu e P R Xi v e - gVx\ Xi Pl e z>* , (2) 

with the 2x2 unitary matrix diagonalising the chargino mass matrix by the right 
(the interactions and notation used can be found in Refs. [9,10], and follow the conven- 
tions of Ref. [11]). By inspection of the Lagrangian it is easily seen that in the decay 
u e — > e~xf, mediated by the first term in Eq. (|2j), the produced electron has negative 
chirality (and thus negative helicity, neglecting the electron mass) and the chargino 
has positive chirality. Since sneutrinos are spinless particles, angular momentum con- 
servation in the v e rest frame implies that the xf must have negative helicity as well, 
as depicted schematically in Fig. [2] (a). For z>* — > e + x{ , the positron has negative 
chirality, and thus positive helicity; therefore, the Xi has positive helicity in the z>* rest 
frame, as shown in Fig. [21(b). 



e v" e e v* 




Figure 1: Feynman diagrams for v* e v e production in e + e annihilation. 



XT e+ Xi 



(a) (b) 
Figure 2: Helicity of the charginos produced in v e and u* decays. 



The decay of Xi to light fermions /, /' and a neutralino Xi is mediated by the 
diagrams in Fig. with similar diagrams for xf decay. We study final states in which 
one of the charginos decays leptonically and the other hadronically. In leptonic decays 
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we restrict ourselves to //' = (ff = z/ At /i + ). The reason is that in decays to e ± 

and a neutrino, the presence of an additional e + e~ pair from sneutrino decays difficults 
the reconstruction of the final state, while in decays to r 1 * 1 the momentum of the 
charged lepton cannot be directly measured. We sum all hadronic decays ff — ud,cs 
iff = ud,cs). 




Figure 3: Feynman diagrams for the decay Xi ff'Xii f° r light fermions //' = 
u e e~,u^~,ud,cs and negligible f L - f R , f L - f R mixing. For / = u e , in diagram (c) 
/l = 1S exchanged. 

We choose a SUSY scenario similar to SPSla but with a heavier sfermion spectrum 
and complex phases <f>i, M . The low-energy parameters most important for our analysis 
are collected in Tabled For vanishing 0i, M , they approximately correspond to rayi = 
250 GeV, m E = mi = niHi = 200 GeV, Ae = —200 GeV at the unification scale, and 
tan/3 = 10. We use SPheno [12] to calculate sparticle masses, mixings and some decay 
widths. Neutralino and chargino masses slightly depend on <f>x, for 0i = M = 
they are m^o = 99 GeV, m~- =178 GeV, m~- = 401 GeV. The relevant branching 
ratios (taking 4>i = (f)^ = 0) are Br(z> e — > e~ xt) — 0.52, Br(x^ — > v^[i~Xi) = 0.10, 
B r (xr ~^ qq'xT) = 0.34, with the same rates for the charge-conjugate processes. 

In the scenario selected the three-body decays Xi ~^ ff'Xi are mediated by off-shell 
intermediate particles. Still, they are dominated by W exchange, diagram (a) in Fig. OH 
due to the large sfermion masses. In Xi ^V/^ Xi the second diagram in importance 
is sneutrino exchange in diagram El (c), but its contribution to the partial width is 13 
times smaller than the one from W exchange. For hadronic decays, diagrams involving 
squarks are even more suppressed, and their contribution has a relative magnitude 
of 2 x 10 -4 with respect to the W diagram. Therefore, many characteristics of this 
scenario are shared with scenarios with a heavier chargino spectrum, in which the two- 
body decay Xi ~^ W~Xi is possible and diagrams (b) and (c) are negligible both for 
leptonic and hadronic final states. 
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Parameter 


Value 


Mi 


102.0 e^ 1 


M 2 


192.0 




377.5 


tan /3 


10 




252.4 




264.5 




571.5 


m d L , m s L 


577.0 



Table 1: Low-energy parameters (at the scale Mz) for the SUSY scenario used. The 
dimensionful parameters are in GeV. 



For values of 0i, M different from 0, 7r, CP is violated in the neutralino and chargino 
sectors, leading to large SUSY contributions to EDMs. If the squark spectrum (which 
does not play any role in our analysis, as seen from the previous discussion) is heavy 
enough, experimental limits on the neutron and Mercury EDMs can be satisfied. On 
the other hand, for the selectron and sneutrino masses under consideration, the exper- 
imental bound on the electron EDM d e severely constrains the allowed region in the 
(0i, <pn) plane. Using the expressions for d e in Ref. [13] it is found that, for each 0i 
between and 2ir, in this scenario there exist two narrow intervals for M , one with 
values 0^ ~ and the other with values M ~ 7r, in which the neutralino and chargino 
contributions to d e partially cancel, resulting in a value compatible with the experi- 
mental limit d e exp = (0.079 ± 0.074) x 10" 26 e cm [14]. For instance, for fa = vr/2 the 
phase M can be M ~ —0.12 or M ~ 3.21. Several representative examples of these 
pairs of phases allowed by EDM constraints are collected in Table[2j We can see that in 
principle it is possible to have any phase 0i, though with a strong correlation with </> M . 
If 0i and M are experimentally found to be non- vanishing, a satisfactory explanation 
will be necessary for this correlation, which apparently would be a "fine tuning" of their 
values [15]. 



3 Generation and reconstruction of the signals 

The matrix element of the resonant processes in Eq. (|T]) are calculated using HELAS 
[16], including all spin correlations and finite width effects. The relevant terms of 
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01 




01 


0M 








7T 





tt/8 


-0.0476 


7tt/8 


-0.0454 


tt/4 


-0.0876 


3tt/4 


-0.0845 


3tt/8 


-0.1136 


5tt/8 


-0.1114 


tt/2 


-0.1218 







Table 2: Examples of approximate phases 0i, 0^ for which the pairs 
(—01, — (f)fj) are allowed by EDM constraints. 



the Lagrangian can be found in Refs. [9, 10]. We assume a CM energy of 800 GeV, 
with electron polarisation P e - = —0.8 and positron polarisation P e + = 0.6. Beam 
polarisation has no effect on the angular distributions and asymmetries studied but 
increases the signal cross section. The luminosity is taken as 534 fb _1 per year [17]. In 
our calculations we take into account the effects of initial state radiation (ISR) [18] and 
beamstrahlung [19,20]. For the design luminosity at 800 GeV we use the parameters 
T = 0.09, N = 1.51 [17]. The actual expressions for ISR and beamstrahlung used in 
our calculation are collected in Ref. [9]. We also include a beam energy spread of 1%. 

We simulate the calorimeter and tracking resolution of the detector by performing 
a Gaussian smearing of the energies of electrons (e), muons (/z) and jets (j), using the 
specifications in Ref. [21], 

ze< io%„ l% _ ^ = om%E ^ A* «e4%, (3) 



E> X 'W I-- 11 ' Ei ^Ei 

where the two terms are added in quadrature and the energies are in GeV. We apply 
kinematical cuts on transverse momenta, px > 10 GeV, and pseudorapidities \rj\ < 2.5, 
the latter corresponding to polar angles 10° < 9 < 170°. We also reject events in 
which the leptons or jets are not isolated, requiring a "lego-plot" separation AR = 
a/ Ai] 2 + A0 2 > 0.4. For the Monte Carlo integration in 10-body phase space we use 
RAMB0 [22]. 

For the precise measurement of angular distributions it is crucial to reconstruct 
accurately the final state. This is more difficult when several particles escape detection, 
as it happens in our case. Since we are not interested in the mass distributions we 
use as input all the sparticle masses involved, which we assume measured in other 
processes [3]. Let us label the electron and positron 4-momenta as p e - , p e + , respectively, 
and the momenta of the "visible" chargino decay products (the ^ for the leptonic 
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decay and the quark-antiquark pair in the hadronic decay) as py+ , Pv- ■ The unknown 
momenta of the "invisible" chargino decay products (the vxi pair in leptonic decays and 
the neutralino in hadronic decays) are pj+, pj- (8 unknowns). From four-momentum 
conservation and the kinematics of the process we have the relations 



where Eqm is the CM energy and pj the missing momentum. These 8 equations are 
at most quadratic in the unknown momenta pi+ , pi- . With a little algebra, they can 
be written in the form of 7 linear plus one quadratic equation, which can be solved 
yielding 2 possible solutions for the unknown momenta (note that in leptonic decays 
only the sum of the neutrino and neutralino momenta can be determined). In order 
to select the correct one we use the additional constraint that either p 2 I+ = m| or 
pj_ = m|o depending on which chargino decays hadronically (this is decided event 
by event depending on the charge of the muon). With pj+, pj- obtained from the 
reconstruction process, the momenta of the two charginos and the two sneutrinos can 
be determined, as well as their respective rest frames. 

We note that ISR, beamstrahlung, particle width effects and detector resolution 
degrade the determination of pi+, pi-, being ISR and beamstrahlung the most trouble- 
some. Detector resolution affects the measurement of charged lepton and jet momenta, 
while ISR and beamstrahlung modify the beam energies and thus reduce the CM en- 
ergy, causing also that the CM and laboratory frames do not coincide. Moreover, the 
last four of Eqs. (J3J only hold for strictly on-shell sneutrinos and charginos. Due to 
these effects, Eqs. (J1J) sometimes do not have a real solution, i.e. the discriminant of 
the quadratic equation mentioned above is negative. In such case, we force the system 
to have a real solution by setting the discriminant to zero, what has the consequence 
that the solutions do not fulfill Eqs. (JH) for the input values m~- , mj> e used but rather 
for other (sometimes very different) ones. 

The determination of the unknown momenta is done as follows. In order to par- 
tially take into account the decrease in the CM energy we replace Eqm in Eq. (|H) by 



E I+ + £ f - 
Pi- + Pi+ 




(Pe 
{Pe 



:+ 



{pv+ +pi+Y 
(pv- +pi-Y 

+ Pv+ +Pi+) 2 

+ pv- +pi-Y 



(4) 



8 



an "effective" CM energy -E^m- We try the reconstruction of the unknown momenta 
for different energies E^ A < E CM in decreasing order, choosing the one which gives 
reconstructed sneutrino, chargino and neutralino masses closest to their true values. 
In case that different effective energies yield equal results for the reconstructed masses, 
we choose the largest one. If the event does not reasonably fit into the kinematics 
assumed for the process, it is discarded. 

For the analysis of some observables we take advantage of the good c tagging ca- 
pability that it is expected at a future linear collider. We use a c tagging efficiency of 
50% and a mistag rate of 0.2% [23]. The identification of c against c, when needed, 
can be indirectly done looking at the charge of the muon produced. In final states with 
H~ , the chargino decaying hadronically is xti thus the tagged jet corresponds to a c 
quark and the other one to a s antiquark. By the same argument, for /i + final states 
the tagged jet corresponds to a c antiquark and the other one to a s quark. 

The backgrounds to the processes in Eq. JH) are rather small. Other SUSY particle 
production processes leading to a final state of e + e~fi ± jj plus large missing energy and 
momentum are for instance e + e~ — > XiXli ~~ *" Xi ^X\i e + e~ — > X1X2 ~* Xi W^XiL 
e + e~ — > X2X34 ~^ X2xfW T } with subsequent decays Z — > e + e~, x\ — > e + e~x1 an d 
hadronic or leptonic Xi 3 decays giving a muon, two jets, a neutrino and a neu- 
tralino. The total cross section of the three processes amounts to 0.1 fb. More im- 
portant is the SM background e + e _ — ► e + e~ n^vqq' (which includes on-shell ZW + W~ 
production) plus its charge conjugate, with a cross section around 4 fb [24]. All these 
backgrounds are expected to be considerably reduced by the signal reconstruction pro- 
cess, which requires that the kinematics of the events is compatible with sneutrino pair 
production. Six fermion production may be further suppressed with a kinematical cut 
requiring that the invariant mass of the e + e~ pair is not consistent with Mz- 

4 Results 

We present our numerical results taking <j>i = M = everywhere except for the study 
of CP asymmetries. We collect in Table El the total cross section for the processes in 
Eq. JH) (including ISR, beamstrahlung and beam spread corrections), the cross section 
after "detector" cuts, its value including also reconstruction cuts, and finally requiring 
one c tag as well. In the latter case, only chargino decays to qq' = cs (qq' = cs) 
contribute in practice. 
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Cross section 



Total 



17.56 



Detector 



11.05 



Reconstruction 



9.99 



c tagging 



2.50 



Table 3: Total cross section (in fb) of the processes in Eq. (|TJ) before and after detector 
and reconstruction cuts, and including also c tagging. 

We give our theoretical predictions for angular distributions (calculated with suf- 
ficiently high Monte Carlo statistics) together with a possible experimental result for 
one year of running (with an integrated luminosity of 534 fb^ 1 ). The latter is generated 
as follows: for each bin i we calculate the expected number of events in one year 
then the "observed" number of events n« in each bin is randomly obtained according 
to a Poisson distribution with mean //j. This procedure is done independently for each 
kinematical distribution studied. In order to be not too optimistic in our results, we 
present samples of possible experimental measurements in which the total number of 
events is approximately 1 a away from the theoretical expectation. 

4.1 Electron angular distributions 

Sneutrino decays are predicted to be isotropic in their rest frame, as corresponds to 
spinless particles. This fact can be experimentally tested in the processes discussed 
here, with the examination of the e~ (e + ) angular distribution in the v e (P*) rest 
frame. We concentrate on the first case, defining /3 e -, as the angle between the e~ 
momentum (in v e rest frame) and an axis orthogonal to the beam direction arbitrarily 
chosen. In Fig. B]we show the dependence of the cross section on this angle. The full 
line corresponds to the theoretical prediction, which slightly deviates from a flat line 
due to detector and reconstruction cuts. The points represent a possible experimental 
result. Despite the small statistical fluctuations, it is clear that the result corresponds 
to a flat distribution. Performing the analysis for three orthogonal axes shows that 
the v e decay is isotropic, what provides a strong indication that sneutrinos are scalars 
and thus that charginos have half-integer spin (as implied by total angular momentum 
conservation). 



10 



400 r 

350 - 
300- 
250 - 
§200- 

u 

150 
100 
50- 



-1 



-0.5 





cos p - 



MC th. - 
Exp. 
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Figure 4: Angular distribution of the e with respect to an axis orthogonal to the beam 
direction, in the v e rest frame. 



4.2 Angular distributions of Xi decay products 

The reconstruction of the chargino and sneutrino momenta allows for the determina- 
tion of the chargino spin directions, which are s_ = P%-, s+ — ~^x + ^ or 
respectively, with P--, the chargino three-momenta in the z> e , z>* rest frames. The 
knowledge of the chargino spin directions then makes it possible a precise study of its 
polarised differential decay widths. The analytical expressions for these quantities are 
rather lengthy [25]; however, the angular distribution of a single decay product in the 
X\ rest frame can be cast in a compact form. Let us define 9j, 9f, 9 as the angles 
between the three-momenta of / = u,u,c, f = n~,d,s and xi m the Xi res t frame, 
respectively, and the Xi spin ( see Fig. EJ). Analogous definitions hold for the xt decay 
products. Integrating all variables except 6j, 9f or 9 , the angular decay distributions 
read 

1 1 + hf cos 9 j 



rHdcos^/ 2 
1 dT(-) l + hfcasOfi 



r(-)dcos9 f , 2 

(5) 



1 dTH l + 4 _) cos0 o 



r(-)rfcos^ 2 
with r^ - ^ = V(xi — > ff'Xi)- The /i factors are called "spin analysing power" of the 
corresponding fermion /, /', Xi, an d are constants between —1 and 1 which depend on 
the type of fermion and the supersymmetric scenario considered. For xt decays the 
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angular distributions dT^ + > /d cos 9 f p are given by similar expressions with constants 
hf, hp, which satisfy hf = —hp hp = —hp, h^ = —h^ if CP is conserved. 




Figure 5: Definition of the polar angles Of, 9p and 9 in the decay of Xi 



The angular distribution of the muon resulting from Xi decay is shown in Fig. 
For cos0 M - ~ —1 the muon is produced opposite to the Xi momentum, then approxi- 
mately parallel to the e + momentum (up to a Lorentz boost). Then, these events are 
suppressed by the requirement of lego-plot separation AR > 0.4. With an accurate 
modelling of the real detector, all the cos 6^- range could be eventually included in a 
fit to experimental data. In our study we restrict ourselves to the ranges where kine- 
matical cuts do not alter the distributions. The fit to the data points (discarding the 
first three bins) gives h^- = —0.270 ± 0.016, in good agreement with the theoretical 
value V = -0.252. 

220 
200 
180 
160 
140 
1 120 

CD 

ca ioo 

80 
60 
40 
20 

°-l -0.5 0.5 1 

cos 9 - 
n 

Figure 6: Angular distribution of the fi~ with respect to the Xi spin, in the Xi res t 
frame. 
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The study of quark angular distributions requires c tagging to distinguish between 
the two quark jets. This reduces the signal by a factor of four, and thus decreases 
the statistical accuracy of the measurements. The angular distribution of the s quark 
is presented in Fig. The suppression around cos# s ~ — 1 is again caused by the 
requirement of lego-plot separation. The fit to the data points gives h s = —0.151 ± 
0.020, to be compared with the real value h s = —0.149. The distribution of the c 
antiquark is presented in Fig. El In addition to the suppresion at cos^g ~ —1, we 
notice a decrease around cos^g ~ 1, which is indirectly caused by the depression at 
cos9 s ~ —1. We thus discard the last seven bins for the fit to the distribution and 
obtain h s = 0.387 ± 0.044, being the real value h 5 = 0.339. 

50 



40 
§30 

pq 

20 



'"-1 -0.5 0.5 1 

cos 8 

s 

Figure 7: Angular distribution of the s quark with respect to the Xi spin, in the Xi 
rest frame. 

As shown by these examples, the "spin analysing power" constants governing the 
angular distributions of Xi are expected to be measured with an accuracy which ranges 
from 6% for h^- to 13% for h s . We note that the samples of "experimental measure- 
ments" generated have been chosen so that the total number of observed events is 
around 1 a away from the expected number, thus the difference between the measured 
and true values of the h constants is due to statistics. The analysis of \t decays 
can be carried out analogously. However, it is interesting to point out that in CP- 
conserving scenarios the decays Xi ff'Xii Xi ~~ > ff'Xi can be- summed, substi- 
tuting cos0 M + |C)g — > — cos^-^s, improving the statistics by a factor ~ \/2. Even in 
CP- violating scenarios this is likely to be a good approximation, because h^- + h^, 
h c + h^, h s + fig are negligible at the tree level, as will be shown in the next subsection, 
and get nonzero values only through higher order corrections. 




13 




Figure 8: Angular distribution of the c antiquark with respect to the Xi spin, in the 
Xi rest frame. 



4.3 CP violation in chargino decays 

The construction of CP-violating observables involves the comparison of Xi an d Xi 
decays. Under CP, the vectors in Fig. El transform to 

s_ -> s + , Pf->-pf , f>f-> -Pf , P^o -> -p^o , (6) 

defined for the xf decay. Then, if we sum Xi an d Xi decays the products 

Q 



i* - s ± ■ Pn± > 

Qc = S± • p c (c) , 
Qs = S± • Ps( s ) 



(7) 



are CP-odd. In hadronic decays an additional product s± • p^o can be measured but 
it equals — (Q c + Q s ). Since these products are even under naive T reversal, the CP- 
violating asymmetries 

N(Q X > 0) - N(Q X < 0) 



A 



x 



(8) 



N(Q X > 0) + N(Q X < 0) ' 

with X = //, c, s and iV standing for the number of events, are also T-even. Hence, 
they need the presence of absorptive CP-conserving phases in the amplitudes in order 
to be nonvanishing. In the processes under consideration these phases are generated 
at the tree level by intermediate particle widths, but they are extremely small, giving 
asymmetries A x ~ 1CT 3 . At one loop level, these asymmetries can result from the 
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interference between a dominant tree-level and a subleading one-loop diagram with a 
CP-conserving phase, and are expected to be very small as well. 

Neglecting for the moment all particle widths (which are kept in all our numerical 
computations), the asymmetries A x vanish at the tree level. Using the fact that in 
this approximation the tree-level partial rates T(xi — > ff'xi) an d r(xl ~* ff'Xi) are 
equal (because partial rate CP asymmetries are also T-even), these asymmetries are 

A X = ^±^. (9) 

Thus, hp- h c + h 5l h s + h s vanish at the tree level if particle widths are neglected, 

and take values O(10~ 3 ) (much smaller than the precision of the measurement of the 
individual constants) when particle widths are included. 

A CP-odd, T-odd asymmetry can be built from the product 

Q12 = a± ■ (p qi xp Q2 ) , (10) 

where pq x and p q2 are the momenta of the antiquark and quark resulting from the 
hadronic Xi decay, respectively. (A triple product s± ■ (px 1 x px 2 ), with X\ = c, c, 
X 2 = s, s distinguished by quark flavour instead of baryon number, is CP-even.) The 
quark and antiquark are distinguished using c tagging and looking at the charge of the 
muon, as explained in section [3 The CP asymmetry 

A _ N(Q 12 > 0) - N(Q 12 < 0) 
12 N(Q 12 > 0) + N(Q 12 < 0) 1 ' 

can be sizeable already at the tree level, and without the need of interference between 
different Feynman diagrams. As discussed in section in the SUSY scenario selected 
the Xi decays are strongly dominated by diagram (a) in Fig. 01 A large triple-product 
CP asymmetry is possible because the polarised decay width for Xi — > "~ * ff'Xi 
contains a term proportional to e^ vpcT st P^- PjPJi m x\ Im O^O^ 1 *, where 0^ R denote 
the left- and right-handed parts of the (complex) WxiXi coupling, respectively, and 
the momenta follow obvious notation. In the Xi res t frame we have P^- = (m--,0), 
Pj = pt, Pfi = pf> and this term reduces to a triple product. 

The theoretical value of A\ 2 as a function of 0i, taking 0^ = 0, is shown in Fig. E3 

(a) , where we observe that the asymmetry can reach values O(0. 1) for large phases 
0i ~ 2. For large M values the asymmetry could be large too, as can be seen in Fig. EJ 

(b) . However, for the range |0^| < 0.12 (modulo tt) allowed by EDM constraints the 
variation of A\ 2 is only between ±0.01. 
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(a) (b) 

Figure 9: Theoretical value of An as a function of 0i, for M = (a) and as a function 
of M , for 0i = (b). 



We present in Fig. HH the value of the asymmetry after signal reconstruction, in- 
cluding the corrections from ISR, beamstrahlung, etc. discussed in section El The 
shaded region represents the statistical error for one year of running with an integrated 
luminosity of 534 fb _1 . The variation of the cross section with 0i is not relevant. In this 
plot we have chosen pairs of phases (0i, M ) allowed by EDM constraints: for each 0i, 
we take a 0^ value, with |0 M | < 0.12, for which the neutralino and chargino contribu- 
tions to the electron EDM cancel. Then, we calculate the asymmetry for this (0i,0 M ) 
pair. (Several representative examples of these allowed pairs of phases are collected 
in Table EJ) The maximum differences in the asymmetry between taking 0^ = and 
taking the M values required by EDM constraints are of 10%, found for 0i ~ 0, ir. 



5 Conclusions 

In this paper we have shown that, in SUSY scenarios where the decays v e — > e~xfi ~~ > 
e + Xi are kinematically allowed, sneutrino pair production provides a copious source of 
100% polarised charginos, in which the kinematics of the process allows to reconstruct 
their four-momenta. The large cross section for the full process e + e~ — > z>*z> e — > 
e + Xi e~Xi ~^ e + y (1 /i"Xi e ~QQ f Xi (plus its charge conjugate) and the low backgrounds 
allow precise measurements of spin-related quantities like angular distributions and 
triple product CP-violating asymmetries. 

The e~ distribution in the v e rest frame shows that the v e decay is isotropic, what 
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Figure 10: Asymmetry A12 as a function of the phase 0i, for M values fulfilling EDM 
constraints, as explained in the text. The shaded region represents the statistical error 
for one year of running. 

is a strong indication that sneutrinos are scalars. Since the spin of electrons is 1/2, 
total angular momentum conservation implies that the charginos have half-integer spin, 
being 1/2 the most obvious possibility. Angular distributions in Xi res t frame allow 
to measure the values of the "spin analysing power" constants h^- , h s , h s controlling 
the angular distributions of its decay products, with a precision between 6% and 13% 
for only one year of running. Since the sums h^- + h^+, h c + h 5 , h s + h s are zero if CP 
is conserved, and are expected to be small in CP-violating scenarios (even including 
radiative corrections), the data from xt decays could eventually be included in the 
analyses, improving the statistical precision. 

Despite the fact that CP-violating phases could be detected in CP-conserving quan- 
tities [26], the direct observation of supersymmetric CP violation is extremely impor- 
tant. We have shown that in xt hadronic decays the CP- violating asymmetry in the 
product s± ■ (pg ± x p q2 ) is very sensitive to the phase of Mi, and can have values up 
to A u = 0.07 for 4>i ~ 2. Such asymmetry could be observed with 2.6 o statistical 
significance in one year of running. It is worth comparing this sensitivity with the ones 
obtained in other processes within the same SUSY scenario. In selectron cascade de- 
cays tL — > ex2 ~^ e Xi/ i+ / u ~) a CP asymmetry in the triple product s- (p^- xp^+) can be 
built, with s the x° s P m [27]. In X1X2 production e + e~ — > X1X2 ~^ XiXi^ + ^~, & — e >£*j 
an analogous asymmetry in the product p e + ■ (pt- x p e +) can be defined [10,28]. The 
values of these asymmetries as a function of 4>i (including ISR, beamstrahlung, beam 
spread and detector effects, as well as backgrounds) are shown in Fig. [TT] (adapted from 
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Refs. [10,27], where two years of running are considered instead of one). Comparing 
Figs.[ini[Il]it is apparent that the observation of CP violating effects is much easier in 
sneutrino cascade decays. For large phases 0i — 2 the statistical significance of the CP 
asymmetry A\ 2 is a factor of two larger than for the asymmetries in selectron cascade 
decays and neutralino pair production. Therefore, provided a good c tagging efficiency 
is achieved, sneutrino cascade decays provide a much more sensitive tool to test CP 
violation in the neutralino sector than neutralino production and decay proceses. 



0.010 - 




Figure 11: Triple product CP- violating asymmetries defined in (a) selectron cascade 
decays [27]; (b) X1X2 production [10], in the same SUSY scenario studied here. The 
shaded regions represent the statistical error for one year of running. 
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